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REPORT No. 1. 

PART 1. 



EXPERIMENTAL ANALYSIS OF INHERENT LONGITUDINAL 
STABILITY FOR A TYPICAL BIPLANE. 

By Jerome C. Hunsaker. 



Article 1. 
INTRODUCTION. 

A model of span 18 inches, representing a typical military tractor 
biplane, was tested in the wind tunnel of the Massachusetts Institute 
of Technology. The lift, drift, and pitching moment were measured 
for a series of angles of incidence corresponding to the maximum 
possible changes of flight attitude. Only the discussion of sym- 
metrical or longitudinal changes is given here. A report on the 
lateral stability of the same model is reserved for a later date. From 
the observed rate of variation of the forces and pitching moment, it 
was possible to calculate the "derivatives" needed in the complete 
theory of longitudinal stability in still air. The damping of the 
pitching oscillation was also determined experimentally. 

The method followed is that of L. Bairstow in his extension of 
Bryan's theory. Notation also follows Bairstow. The value of 
Routh's discriminant, which Bryan has shown to be a measure of 
dynamical longitudinal stability, has been calculated for six speeds, 
ranging from the maximum to the minimum possible speeds for the 
aeroplane type selected. The principal point of interest brought 
out in this connection is that stability falls off rapidly as speed 
decreases or angle of attack increases, and that while this aeroplane 
appears to be very stable at high speeds, it is frankly unstable at 
speeds below 47 miles per hour. 

This instability at low speeds takes the form of an oscillation in 
pitch combined with changing in forward speed and a rising and 
sinking of the whole aeroplane, which, therefore, follows an undulatory 
flight path. The period of the undulation is about 12 seconds, and 
the amplitude doubles itself in less than 20 seconds. Obviously, the 
pilot can not safely abandon his controls at slow speed. 

The importance of this demonstrated instability at low speeds 
should be appreciated in view of recent accidents with military 
aeroplanes when operated at slow speeds. 

The entire investigation of inherent longitudinal stability was pre- 
liminary to the discussion of the effect of wind gusts. Naturally, it 
was first necessary to find a stable aeroplane and to obtain some idea 
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of the "range" of stability. It now appears that a typical aeroplane 
is inherently stable in the sense defined at high speeds only. The 
effect of gusts on the uncontrolled aeroplane will, therefore, be 
investigated only for the high-speed condition. At low speeds the 
aeroplane can not be left to itself in still air. Consequently, a dis- 
cussion of its certain destruction if abandoned in gusty air appears 
unprofitable. 

Article 2. 
MODEL AND PROTOTYPE. 

The type of aeroplane selected is a high-speed military biplane 
tractor known as Curtiss JN2. Shop plans of this aeroplane were 
kindly furnished by the Curtiss Aeroplane Co., Buffalo, N. Y., to 
whom acknowledgment must be made for much valuable assistance, 
including the experimental determination of moments of inertia etc 
by Dr. A. F. Zahm of that company. 

The principal dimensions of the aeroplane were assumed as follows: 



Weight full load pounds. . 1, 800 

Brake horsepower horsepower. . 110 

Maximum speed for calculations miles per hour 79 

Minimum speed for calculations do. . 43 7 

Total wing area (including ailerons) square feet! ! 384* 0 

Area fixed tail do . . 23 0 

. Area horizontal rudder '.'.'.'.'.'.'.do'." 1.9 0 

Area vertical rudder . . . \ \ do 7 8 

Span of wings " ".W" .".feet! ! 36.' 0 

Chord of wings. do.... 5.3 

Gap between wings do 5 3 

Length of body '.'.'.60'.'.'.'. 26.' 0 



The model was made geometrically similar to its prototype and 
one twenty-fourth scale. The general features are shown in the 
drawings of the model. (Figs. 1 a, 6, c.) The model was an exact 
copy of the aeroplane except for the propeller and wing wirin°- 
which features were omitted. Also wing struts were made rouiKi 
instead of "stream-line" in section. Since it is well known that 
the resistance of a series of similar aeroplanes varies somewhat less 
rapidly than the square of the speed and square of a linear dimen- 
sion, due to skin friction, it is believed that the prediction of the 
resistance of the full size aeroplane from the observed model resistance 
will still be a fair estimate in spite of omissions on the model. 

For simplicity, the model was made with the trailing ailerons or 
wing flaps integral with the wings. This somewhat increases the 
effective supporting area. Also the fixed tail and elevator were 
made in one, corresponding to the elevator held fast in its neutral 
position. These points are made clear on the drawings of the model. 

Article 3. 

GENERAL WIND TUNNEL PROCEDURE. 

The model was tested in the 4-foot wind tunnel at a velocity of 30 
miles per hour. The wind tunnel and aerodynamical balance are 
duplicates of the installation of the National Physical Laboratory, Ted- 
dington, England, and reference should be made to the Technical 
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Report of the Advisory Committee for Aeronautics, London, 1912-13, 
for detail description and methods of operation. 

In general, it may be stated that the wind tunnel provides a wind 
constant in velocity within 1 per cent, which velocity is further con- 
stant across the working cross section of the tunnel within l£ per 
cent. Velocity is measured by a suction plate calibrated against a 
standard Pitot tube with a precision of one-half per cent. The 
model is mounted on the balance in various attitudes of pitch or yaw, 
and in such positions are measured the three forces and three couples 
produced by the wind along and about three mutually perpendicular 
axes in space. From a knowledge of the variation of these forces and 
coimles with change of attitude, the so-called " resistance derivatives " 
of Bryan's 1 theory of dynamical stability may be computed. 

The theory of stability also requires the determination of the damp- 
ing of oscillations about the center of gravity of the aeroplane. A 
special oscillating apparatus was built for these tests which will be 
described below. By oscillating the model in the wind and observ ing 
the decrement of amplitude with time, it was possible to estimate the 
"rotary derivatives. 

Article 4. 
LONGITUDINAL TESTS. 

The model was mounted on the balance with its wings in a vertical 
plane by means of a vertical rod driven into the body at the point 
shown on figure 1&. By swinging the model about the vertical axis 
passing through the spindle, the angle of wind to the wing chord was 
varied from + 20° to — 8°. At each attitude the force across the wind 
or "lift," force down wind or "Drift," and the pitching moment 
about the spindle were measured. The signs were taken so that 
an actual lift, actual head resistance, and a stalling moment are posi- 
tive. The wind velocity was 30 miles per hour of standard dry air 
at 15° C. and 776 mm. Hg. The experimental points are shown on 
figure 2, w'here forces are in pounds and moments in inch-pounds. 
The precision of measurement is within 1 per cent. 

For a given attitude, the res ultant fo rce on the model in pounds 
at 30 miles per hour is i? = VZ, 2 + Z> 2 . This resultant makes an 

angle with the wind direction given by a = tan~ 1 ^- The force R 

is observed to have a pitching moment M about* the spindle axis. 
It may then be assumed to be situated so that the perpendicular 

from this axis to R is given by z = ^r- The vector R is thus deter- 
mined in magnitude, direction, and line of application. The resultant 
force vectors R are shown on figure 1& to a scale 1 hich equals 0.2 
pound. The vector R is purely an algebraic substitution for the 
complicated system of forces and couples acting on the aeroplane. 
The vectors are drawn relative to the aeroplane. 

The center of gravity was assumed to lie as shown near the inter- 
section of the propeller axis with the resultant force vector for 4°. 
At this attitude, then, the pitching moment should be nearly zero. 

1 G. H. Bryan, Stability in Aviation. 
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The c. g. location determined for the actual aeroplane after exten- 
sive trial flights is almost identical. 

It is seen that for angles smaller than 4°, R passes forward of the 




c. g. and for angles greater than 4°, it passes to the rear. The 
aeroplane is longitudinally stable in a static sense. It will be shown 
below that it is not always dynamically stable. 
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Article 5. 

PERFORMANCE CURVES. 

The lift and head resistance or "drift" of the full scale aeroplane 
were assumed to be approximately given by the relation: 

Force on mode l ( 30 V 
Force on aeroplane \24 V/ 



figure: 3. 
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when V is the flying speed of the aeroplane in miles per hour. The 
above relation holds, of course, only for the same attitude of model 
and aeroplane. The weight of the aeroplane, 1,800 pounds, must 
equal the lift in flight. Hence : 
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A series of speeds V was computed for a series of attitudes of the 
aeroplane, and the aeroplane drift at each attitude was then com- 
puted from: 

2 / y\2 

D full size = Z> model x 24 x( ^ ) 

In figure 3 are given curves of drift, effective horsepower required, 
angle of w r ing chord to wind and ratio weight to drift plotted on V as 
abscissae. For our calculations a maximum speed of 79 miles and a 
minimum of 43.7 miles were selected corresponding to angles of wing 
chord to wind of 1° and 15.5°, respectively. 

The curve of E.H.P. on figure 3, indicates that 87 propeller horse- 
power is necessary for a speed of 79 miles. If the propeller has an 
efficiency of 80 per cent, the motor must develop at least 110 brake 
horsepower. The original designs contemplated as maximum speed 
of about 80 miles per hour for a 120 brake horsepower motor, which 
appears very reasonable. As actually built this type was given a 
rated 90 horsepower motor. Assuming 70 E.H.P. delivered to the 
propeller a speed of 73 miles per hour is indicated by our curves. It 
is reported that the speed of this aeroplane was actually 73 miles per 
hour. 

Article 6. 

CHOICE OF AXES— NOTATION— UNITS. 

Axes for reference are assumed fixed in the aeroplane and moving 
with it in space. The origin is at the center of gravity. < For steady 
horizontal flight at a given attitude the axis of Z is vertical, the axis 
of X Horizontal and directed to the rear in the plane of symmetry, 
and the axis of Y is horizontal and directed toward the left-wing tip. 
Forces along these axes are denoted by X, Y, Z and are expressed in 
pounds per unit mass. Moments are Z, M, N and are given in pounds- 
feet per unit mass. 1 

Angles of roll, pitch and yaw from the normal flying attitude are 
denoted fry f, 6 and ty. Angular velocities of roll, pitch and yaw are 
p, q, r in radius per second. The signs of moments, angles and angular 
velocity are positive considered in the directions XY, YZ or ZX. 

Moments of inertia referred to axes X, Y, Z are denoted by 
mK 2 Al mK 2 B, mK 2 c where m is the mass of the aeroplane and Ka, 
Kb, K c corresponding radii of gyration. 

Article 7. 
EQUILIBRIUM CONDITIONS. 

In normal horizontal flight in still air a state of equilibrium is 
assumed such that the power available maintains the aeroplane at 
such a speed that the weight is just sustained. Since the ltft of an 
aeroplane whig is also a function of its attitude or angle of attack, 
it is further assumed that the attitude is proper for the speed. In 

i Unit mass is the slug equal to 32.17 pounds weight. 
H.",302°— S. Doc. 268, 64-1 3 
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normal horizontal flight the axis of X is parallel to the apparent wind 
direction and is hence horizontal. Let 8 be the angle ox pitch of the 
aeroplane away from its normal attitude. Then normally 6 is zero. 
Likewise if the aeroplane is in equilibrium in its flight, the angular 
velocity of pitch is zero and also the pitching moment, l/ 0 . 

At high speed, for example 79 miles, the axis of X is horizontal and 
makes an angle of 1° with the wing chord. At low speed, new axes 
are chosen such that the axis of X is still horizontal but makes an 
angle of 15.5° with the whig chord. The axes are fixed by the 
equilibrium conditions for flight and differ for each normal flying 
attitude. Oscillations about the normal flight path when the motion 
is disturbed are referred to the above defined axes which are assumed 
fixed in the aeroplane and moving with it in space. 

The pitching moment curve observed for the model shows zero 
moment for an angle of wing chord of 4.5° and a diving moment at 
larger angles. For slow flight, it is assumed that the pilot by proper 
setting of his horizontal rudder impresses an equal stalling moment 
on the machine so that the net pitching moment is zero. The effect 
is to move the pitching moment curve parallel to itself by the alge- 
braic addition of a stalling moment so that its ordinate has zero value 
for the desired flight attitude. 

Article 8. 

TRANSFORMATION OF AXES. 

It is convenient to measure in the wind tunnel the lift and drift 
about axes always vertical and horizontal in space. For the oscilla- 
tions of the aeroplane it is convenient to consider the forces referred 
to axes fixed in the aeroplane as described above. The transforma- 
tion is effected in the usual way by means of the formula: 

m Z' = L cos 0 + Z? sin 0, 
m X f = D cos Q — L sin 0, 

where 0 is the angle of pitch of the aeroplane away from its normal 
attitude, considered positive for stalling angles. Here L and D are 
lift and drift on the model in pounds, and m X' and m Z' correspond- 
ing forces in pounds along the axes X and Z. The model forces 
Z' , X' are converted to Z, X, full size, by multiplying by the 
square of the speed and linear dimension ratios. The following tables 
carry out the required transformation. 

The pitching moment if is independent of the longitudinal shift of 
axes and varies only as the square of the speed. Curves of X, Z and M 
for the different flight attitudes are plotted on figures 4, 5, 6, 7, 8, 
and 9. The transformation of the moment about the spindle to the 
corresponding moment about the c. g. of the full-size aeroplane is 
given below. 
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{=1°, V=79 miles, m=55.9 slugs. 
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1=15.5°, F=43.7 miles. 
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CONVERSION OF PITCHING MOMENTS. 

mM 8 = moment about spindle in inch pounds on model. 
mM) cg == moment about c. g. in inch pounds on model. 
6=3.04 inches, c. g. forward of spindle. 
a=0.10 inches, c. g. above spindle. 
Axis of X 3.5° to wing chord. 

J/=pitching moment about c. g. full size, full speed, in pounds feet per unit mass. 

m M C g=mM s — mZ'b -?nX / a. 

i=angle of wing chord to wind, degrees. 

6=angle of axis of A'' to wind, degrees. 
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Article 9. 

RESISTANCE DERIVATIVES, LONGITUDINAL. 



Notation follows Bairstow, 1 to whose paper reference should be 
made for the detailed discussion of " derivatives. " In the theory of 
small oscillations, the aerodynamic forces X 0 , Z 01 and pitching 
moment, M 0 , are eliminated by the conditions of equilibrium. In 
disturbed motion, disturbances in normal flying speed and attitude 
cause changes in the quantities, X, Z, and M. 

Let V be the normal flying speed and u, w and q small changes in 
horizontal and vertical velocity components and angular velocity of 
pitch. If the disturbance be small, u, w and g are small with respect 
to U. For example, the function 

X=f(U+u, % q) 

may be expanded into the approximate form 

X = X 0 + uX u 4- wX w + q W q , 

a linear function of the small quantities u, w, q. The coefficients 
X u , X WJ X q are the so-called resistance derivatives of the theory of 



1 Technical Report of the Advisory Committee for Aeronautics, London, 1912-13. 
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small oscillations, and physically represent the slope of a curve of X 
on a base u, w, or q. 

Similarly 

Z = Z 0 + uZ u + wZ w + qZ q 
M= Af 0 + uM u + wM w + qM q 

From the conditions of equilibrium, X 0 is balanced by the pro- 
peller thrust, Z 0 by the pull of gravity or Z 0 =g, and M 0 = o. 
Also, Bairstow has shown that X^ and Z q may be neglected. 
e X u is the rate of change of X with change in forward speed. But 
since X is a function of forward speed squared we may write: 



and 



u AU U 



S. Doc. 208, G4-1. 
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These coefficients may be obtained directly by calculation since 

Drift ^ a ^ Z 0 = g. For example, at 79 miles per hour, U= 

— 115.5 feet per second and Z 0 = 32.2. Then 

„ 2X32.2 _ 7 
-115.5 " ' o0 ' 

Also at 15?5, Z7= —63.8 feet per second and 

_ 2X10 
^ tf ""^63^""- 276 

The derivatives X wy Z w , 2f w represent the effect of a vertical 
component of velocity. From the well-known method of velocity 
composition, the vertical velocity w acts with the horizontal velocity 
U to cause the apparent wind to have an inclination to the horizontal 

Df tan -1 -jj- This inclination is given to the model in the wind 

tunnel, and X, Z, and M measured for various pitch angles. 

an An ( 

But Ad = tan -1 -^=57.3 when Ad is a small angle in degrees. 

. v AX 57.3 AX 
' ' w ~ U ' Ad 

AX 

-jTjr is the slope of a curve of X on pitch angle as base. For example, 
from figure 4, = — ^ and 

Similar formulas are used to compute Z w and M w . It may be noted 
that the method assumes that for small oscillations, hence small 
changes 6, the tangent may be substituted for the actual curve. 
The Emit of validity is obviously the range of pitch angle over which 
the tangent to the curve is not greatly changed. This range is usually 
about 4 to 8 decrees. 

The values of the resistance derivatives calculated in this manner 
will be found tabulated later. 

Article 10. 

DAMPING. 

The damping of pitching about the c. g. is represented by the rotary 

derivative M q . For an angular velocity ^ = g, a damping moment 

qMg is exerted on the aeroplane. 

To measure this aerodynamic dampin & , the special oscillating appa- 
ratus was designed which is shown oy the photograph of figure 10. 
The model is mounted on a massive bracket which pivots about the 
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two points shown. Fore-and-aft arms carry counterweights which 
are adjusted to give a reasonable natural period. The spiral springs 
bear in notches on the arms by means of knife-edged shackles. The 
springs insure that the motion shall be oscillatory. The assumed 
c. g. location of the aeroplane model is arranged to be on the axis of 
rotation. The actual center of gravity of the apparatus is not 
considered ^ > 

Friction is kept small by careful design of the steel pilots, which 
are hardened steel points bearing in tool steel cones. The spring 
knife edges are glass hard. It was found that a convenient period 
is about one-half second. In still air the apparatus will rock more 
than 300 times before the amptitude is diminished by friction to 
one-ninth of the initial displacement. m 

The moment of inertia of the entire oscillating mass was calcu- 
lated and then checked by an independent experimental determina- 
tion. 



Let: 



7= moment of inertia of all oscillating parts in slug foot- 
units. 

ra' = mass of all oscillating parts in slugs. 

M 0 = moment of air forces on model at rest. 

3f, = moment of springs at rest. 

Ed = additional moment of springs when deflected. 

c = c. g. of entire apparatus above pivot, feet. 

0 = angle of pitch from normal attitude in radians. 

Mo -^ = damping moment due to friction. 

p w ^ = damping moment due to wind on apparatus. 

M m J^= damping moment due to wind on model. 
cm' 6 = static moment due to gravity. 



The equation of motion then is: 



But 2f 0 = M„ by the initial condition of equilibrium. Let 

p-M-aH-M then t^r +#£ + OC-«"0* : « 

The solution of this equation is well known to be: 



e - C e ™ eoJ<y ( JT- cm' ) \ - jj 2 + oc 

where 0 and oc are arbitrary constants. If time be counted when 
the amplitude of swing is a maximum then cos{ — } = 1, and QjK 
the initial displacement. Also if the number of beats be counted by 
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observing the times for succeeding maxima, a plot of amplitude oil 
time wilfhave for its equation the simple form: 



-Ml 

e = e 0 e a 



The coefficient p is the logarithmic decrement of the oscillation 
and must be numerically positive to insure that the oscillation dies 
out with time. ■ . , . u 

The apparatus was fitted with a small reflecting prism by wtnch a 
pencil of light was deflected toward a ground glass plate set m the 
roof of the tunnel. Nine lines spaced 0.2 inch were ruled on this 
plate With the model at rest the beam of light was brought to a 
sharo focus on the line marked zero. By means of a trigger the 
observer started an oscillation of the model, and the spot of light was 
observed to oscillate across the scale. The time ,t was observed in 
which an oscillation was damped from an amplitude of 9 to an ampli- 
tude of 1, for example. 

Then: log^= *V=%A and knowing / and t, p is calculated. 
6 2,1 

Preliminary tests showed that the same value of ix was obtained 
whether the timing stopped at 0 = 5, 4, 3, 2, or 1. 

Oscillation tests were made at five wind velocities varying Irom 
5 to 35 miles per hour. The coefficient ii appeared to vary approxi- 
mated as the first power of the velocity. t 

Similar tests were made with the model for no wind to determine 
u„ which may be said to be due almost whoUy to friction and verv 
slightly to the damping of apparatus and model moving through 

Likewise ,x w was obtained by oscillating the apparatus without 
model in winds from 5 to 35 miles per hour. 

The coefficient Mm has the dimensions 1 pYV .where p is density of 
air, 1 a linear dimension, and Fthe velocity of the wind To convert 
u to M„ for the full-size machine at full speed, multiply by the fourth 
power of 24, the scale, and by the ratio of full speed to model speed. 

The numerical results of tests of the pitching oscillation follow. 
Note that the damping of the pitching falls off for low speeds, ihis 
contributes to the difficulty of providing sufficient stability at low 

SP In d the tables following, the number of beats, n, is recorded as a 
general check and is not used. Recorded values of n and t are the 
means of three or five separate observations^ ^ 

' Bairstow, loc. cit., p. 176. 
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PITCHING OSCILLATION TESTS. 

I model and apparatus =0. 04195 
I apparatus = . 0368 

Apparatus. 

Wind velocity, miles. ... 0 14. 7 21. 4 35 

n beats counted 350 253 210 186 

/ seconds 168 120 100 90 

fj. 00096 .00135 .00162 .00180 

Hv> (less zero) 0 . 00039 . 00066 . 00084{ U ^ 

Apparatus and model with wing chord 1° to wind. 

Fniiles 0 9.5 14.7 21.3 25 < 30 37.3 

n beats 300 90 56 40 35 32 27 

t seconds 160 45 28. 5 20 17. 5 16 13. 5 

fx gross 00115 . 00410 . 00646 . 0092 . 0105 . 0115 . 0137 

Mo friction 00096 .00096 .00096 .0010 .0010 .0010 .0010 

ti w apparatus 0 . 00035 . 00040 . 0006 . 0007 . 0009 . 0011 

Mm net 00019 . 00284 . 0054 . 0076 . 0088 . 0096 . 0117 

But pi m = —m M q when reduced to full size and 79 miles per hour 
and mass of 55.9 slugs. 



.'.M q = -.009GX(24) 4 X (79/30) X 1/55.9= -150.0 

or for 

£7= -114 foot-seconds, M q =1.32 U 

Apparatus and model with wing chord 15.5° to wind. 

V. 9.1 14.7 21.4 25 30 

n 75 50 35 30 25 

t 38.5 25.0 17.5 15 13 

/x^ross 0048 .0074 .0105 .0123 .0142 

jx TO net 0035 .0060 .0089 .0106 .0123 



37.5 
19 
9 

.0205 
.0184 



of 



Ma 



.0123 X (24) 4 X (43. 7/30) X 1/55.9 = 



106 



M Q =1.66 U where Z7 is —64 foot-seconds, or 43.7 miles. 



The computed values of p. m , the model damping coefficient, are 
plotted on figure 11. It appears that /% is approximately a linear 
function of the velocity, as would be expected, and the conversion 
to full scale, full speed, is made as indicated above. 

The damping coefficient is not greatly different for different atti- 
tudes, and the following values are obtained by interpolation: 



Angle of 
wing chord 
to wind. 

+1° 
7° 
10° 
12° 
14° 
15. 5° 



V. 
79.0 
51.8 
47.0 
45. 2 
44.2 
43.7 



U. 
-115.5 

- 75.8 

- 68.8 

- 66.2 

- 64.8 
-64.0 



Mq. 

L 30U=-150 
L 49 U=-113 
1. 55 U=-108 
1. 59 U=-106 
1. 63 U=-106 
1. 66 17= -106 
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Article 11. 
RADIUS OF GYRATION. 

For the radii of gyration of the fully loaded aeroplane we are in- 
debted to Dr. A. F. Zahm. The actual aeroplane, complete with 
gasoline, water, pilot, passenger, and other weights in place, was 
suspended from a beam by a chain. The center of gravity was first 
located by an inclining method. The machine was then made to 
oscillate in pitch about the point of attachment of the upper end of 
the chain. Light guys were run to tail and wing tips to insure that 
the chain and aeroplane moved as a rigid body. 

Let the distance from center of gravity to point of suspension be 
denoted by h, p the natural period of oscillation in seconds, K B the 
radius of gyration in feet about the Y axis or axis of pitch, then 

| : "-(£)*-» 

By observation Ji = 12.2 feet, p = 60/14 seconds. 

Z* 2 = 34, K B = 5.S feet. 1 

Article 12. 

ROUTH'S DISCRIMINANT. 

Bryan 2 has shown that the character of the longitudinal motion 
of an aeroplane may be investigated with reference to the roots of 
a biquadratic equation of the form: 

A\*+B\* + C\ 2 + D\ + E= 0 

The equations of motion may be considered of the form Q=Ke u 
where K is some constant. For stability the quantity X must be 
negative if real, or have its real part negative if complex, in order 
that the amplitude of the motion will diminish with time. 

The condition that the real roots and real parts of imaginary roots 
of a biquadratic equation with constant coefficients shall oe negative 
is that the coefficients A, B, C, D, E shall each be positive as well 
as the quantity BCD — AD 2 — B 2 E. The latter is commonly known 
as Routh's 3 discriminant. 

The constant coefficients A, B, C, D, E, are functions of the con- 
stants of the aeroplane at the normal flying attitude, i. e., the follow- 
ing: X Uf X W1 X q , Z u , Z w , Z q , M u , M w , M q , U, and K B 2 . These are 
resistance and rotary derivatives, velocity, and radius of gyration. 
For a given attitude and for small oscillations about that attitude, 
it is considered that these quantities are constant. For simplicity 
it is here assumed that normal flight takes place in a horizontal plane 
and the inclination of the flight path and consequent components of 
gravity in the axes of X andZ are eliminated. Also X q and Z q are 



1 It is of interest to note that the radius of gyration for rolling was estimated to be 6.2 feet. 

• Stability in Aviation. 

• Advanced Rigid Dynamics, E. J. Routh. 
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neglected as unimportant and M u is zero by the conditions of equi- 
librium. For the computation of Routh's discriminant we require to 
know, then, only those quantities which have been so far determined, 
and which are assembled below for the different cases investigated. 

Formulae for the coefficients A, B, C, D, E are given by Bairstow 
and are used here, but making 9, X q , Z q , and zero. They are 
copied in simplified form for reference. 

A=K B 2 

B=-(M q + X u I& + Z w K B >) 
+ X u M q +K B * 



z w , u 

M w , M a 



D = 

E=-gM w Z, 



Z U i Z w , TJ 



Article 13. 



BAIRSTOW'S APPROXIMATE SOLUTION. 

From consideration of the usual relative numerical values of the 
coefficients of the biquadratic, Bairstow has shown that the equation 
may be factored to a first approximation and put into the following 
form: 



(x 2 +5/ix+ c/a) (x 2 +[z>/c-^f] X- 



in which the first factor represents a very short oscillation, which 
in most aeroplanes rapidly dies out and is of no importance. The 
second factor represents a relatively long oscillation involving an 
undulatory flight path with changes in pitch, forward speed, and 
altitude. The long oscillations should diminish in amplitude with 
time, in which case the motion is stable and the aeroplane will return 
to its original normal flight attitude if temporarily deviated there- 
from by accidental cause. The motion is unstable if the long oscilla- 
tion increases in amplitude with time. It will be shown that the 
aeroplane under investigation is stable at high speeds and unstable 
at very low speeds. It is believed that this is true of all aeroplanes. 

Case I. 

i=incidence, wing chord to wind +1°. 

Velocity F=79 miles. Z7= — 115.5 foot-seconds. 
?n=55.9 slugs, ^=34. 



X U -.12S 
Z u - .557 

A=+ 34 
£=+289 
C=+834 
£=+115 
£=+ 31. 



Z w -3.95 



3^+1.74 

M q -150 



BCD-AB*-B 2 E=+1SX10 6 stable. 
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Short oscillation: X 2 +8.5X-f- 24.5=0 
X=-4.25±2.54i 

p=period=^54 =2 -^ seconds. 

/=time to damp 50 per cent=^||=.16 second. 

Long oscillation: X 2 +.125X-h 0374=0 

X=-.063±.183i 
p=34.3 seconds, t =10.8 seconds. 

The short oscillations are unimportant. The long oscillations are 
easy and strongly damped. The aeroplane should be very steady 
at this speed. 

Case II. 

i=7°. F=5 1.8 miles, U= — 75.9 foot-seconds. 

X u -.\2l ^+.113 Jf„+2.45 
Z„-.849 £^-2.26 M q -l\Z 



A = + 34.0 
£=+194.0 
C= +467.0 
D=+ 64.3 
E=+ 67.0 



BCD-AD 2 -B 2 E=+?>2X10 5 stable. 



Short oscillation: X 2 +5.7X+ 15.9=0 
X=-2.85±2.33i 
p=2.7 seconds 

t =.24 second to damp £0 per cent. 

Long oscillation: X 2 +.078X+. 143=0 
X=-.039±.377t 
p=16.7 seconds 

t =17.7 seconds to damp 50 per cent. 

The period is shorter than at high speed and the damping less 0 
The aeroplane should therefore be less comfortable. 

Case III. 

t=10°, F=47 miles, U= -68.8 foot-seconds. 

X U -.VS\ Z„-.936 3/^+2.50 
-X^-,075 £*HM6 M q -108 



A= + 34 
£=+165 
C=+355 
Z>= + 42.5 
E=+ 75.3 



BCD-(AD 2 +B 2 E)=3.$X10 5 stable. 



Short oscillation: X 2 + 4.85X+10.44= 
X= -2.42+2.12 



p=2.96 seconds. 

t = .28 second to damp 50 per cent. 

Lone: oscillation: X 2 + .021X+.212=0 
X=-.011±.460i 

p=13.71 seconds. 

t =62.7 seconds to damp 50 per cent. 

This oscillation is rapid and but slightly damped, and would 
probably be uncomfortaole. The stability is slight and wind gusts 
or external disturbances, if recurrent, might cause trouble. 
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Case IV. 

t=12°, F=45.2 miles Z7=— 66. 2 foot-seconds. 



x u -.m 

X„-.236 



Z u -.972 
Z U i— .736 



Jf g -106 



A = + 34 1 
£=+137.5 
C = +243 
D= + 17.4 
E=+ 67.2 



BCD-AD 2 -B 2 E=-7XW Unstable. 



Short oscillation: X 2 + 4.04X+7.14=0 
X= -2.02+1. 75i 



/;=3.59 seconds. 

/ = .342 second to damp 50 per cent. 

Long oscillation: X 2 - .985X=.276=0 
\ == + .043 ±.524i 

p=12.0 seconds. 

t =16.0 seconds to double amplitude. 

The machine is frankly unstable and the i^ilot dare not release his 
elevator control. 

Case V. 

i=14°. F=44.2 miles, Z7= -64.8 foot-seconds. 



X u -.223 
A' -.132 



Z u -.993 
Z in — .553 



M q -106 



A = + 34 
B =+134 
C=+213 
Z) = + 28 
E=+ 63.6 



BCD -AD 2 - B 2 E= - 3 . 7 X 10 5 Unst ab le . 



Case VI. 

i=15.5° F=43.7 miles, Z7= -63.8 foot-seconds. 



X u -.276 
A\, ; -.292 



Z„-1.01 
Z w — .673 



if^+2.02 

M q -m 



A = + 34 
£ =+138 
C=+226 
Z)=+ 24.2 
E=+ 65.7 



BCD-AD 2 -B 2 E=-bXW Unstable. 



Short oscillation: X 2 + 4.06X+6.65=0 
X=-2.03± 1.59t 



p=3.95 seconds, period. 

t = .34 seconds to damp 50 per cent, 



Long oscillation: X 2 + .071X+.291=0 
X = + .0358zb.541i 



25302°— S. Doc. 268, 64-1- 
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Real part of X is here positive, indicating an oscillation increasing 
with time. 

7>=-^j=11.6 seconds. 

t == ~(^j^j := 19.3 seconds to double amplitude. 

The motion is both rapid in period and rapidly increasing in am- 
plitude. Indeed the amplitude is doubled in two swings. This aero- 
plane, if left to itself, would be highly unstable. 

Article 14. 

VARIATION OF LONGITUDINAL STABILITY WITH SPEED. 

Preliminary to consideration of the action of gusts on an inherently 
stable aeroplane, it was desired to analyze the motion in still air of a 
machine which could be called inherently stable longitudinally. It 
has been found above that a typical aeroplane becomes less stable 
at low speeds until real instability results. This result is somewhat 
unexpected in view of the curves of pitching moment M, which in- 
dicated static stability at all possible attitudes up to and including 
horizontal flight at + 15°. 5. In other words, M w is positive for aft 
cases. The instability comes about on account of the rapid rate of 
increase of drift at large angles causing X w to change sign, and on 
account of the less rapid rate of increase of lift, causing Z ?<7 to be- 
come small at high angles of pitch. Furthermore, M q diminishes at 
low speed. 

From the speed power curves on figure 3, it appears that for angles 
greater than 10°, we are on the part of the power curve which re- 
quires more power to go slower, u region of reversed controls." This 
region is now found to be dynamically unstable so that controlled 
flight only is possible here. But with reversed controls this is 
doubly dangerous. 

The frequency of accidents at low speeds, following the recent 
demand for a wide speed range, confirms this impression of the 
danger of low speeds when approaching a critical angle and speed. 
The critical angle for instability is clearly an angle less than the pos- 
sible maximum for flight. 

A fair measure of tne relative stability at various speeds may be 
had by noting the following tabulation of the values of Routh's 
discriminant, denoted by R: 



Velocity 


Wind chord 




in 


to 


B. 


miles. 


wind. 




79.0 


1° 


+180 X10 5 ] 


51.8 


7° 


+ 32 X Unstable. 


47.0 


10° 


+ 3.8X10 5 J 


45.2 


12° 


- 7 X10 5 ) 


44.2 


14° 


- 3. 7 X10 5 [Unstable. 


43.7 


15. 5° 


- 5 X10 5 J 



The table is reproduced graphically on figure 12. 
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A similar investigation for lateral stability fails to show any marked 
change with speed, as would be expected since speed depends on 
pitch angle and the factors which make or unmake lateral stability 
are but slightly affected by angle of pitch. 
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VARIATION WITH ATTITUDC. 



REPORT No. 1. 

PART 2. 



THEORY OF AN AEROPLANE ENCOUNTERING GUSTS. 

By Edwin Bid well Wilson. 



Article 1. 
INTRODUCTION. 

Z Sfe P w pages ™ shall begin at the start «*h some 

fer'Iitiass'oHh'e Lropiane); then the dynanucal equations of 

motionare vm***-"-** (la) 

dw'/dt + v'p'-u'q'=Z', W 



dhJdt-r% + q% = mL', 
dhJdt-p'h+r'K^mM', 

dh 3 /dt-q%+P%= mN '> 

where m is the mass and 

h^p'A-q'F-r'E, 

Ti^q'B-r'D-p'F, 

l^r'C-p'E-q'D, 



(2a) 
(26) 
(2c) 

(3a) 

m 

(3c) 
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are the components of angular momentum, — the quantities A, B, G 
being the moments and D, E, F the products of inertia relative to 
the moving axes fixed in the body. 

The symmetric aeroplane will alone be considered here; 

D-F=0. (4) 

If the machine is in uniform horizontal flight, all the forces, 
moments, linear velocities and angular velocities except u' vanish, 
and u' = TJ, a negative quantity in magnitude equal to the uniform 
velocity. (The precise backward direction of the z-axis is that 
which is horizontal in uniform flight, and hence by this definition the 
direction of this axis, and of the z-axis, varies in the aeroplane with 
the speed.) 

If the motion is slightly disturbed, the velocities take the values 

u' = U+ u, v' = v, w f = w, p r = pj q f = q, v' = r, (5) 

where u> v y w, p, q, r are small. The products of these small quan- 
tities are neglected, as in all discussions of small oscillations, and the 
equations take the form 

du/dt = X' y dv/dt + Ur= Y', dw/dt- Uq = Z', (6) 

Adp/dt - Edr/dt = mU , Bdq/dt = m W , Cdr/dt - Edp/dt = mW . (7) 

In uniform motion the forces and moments all vanish. For the 
disturbed motion they are small and may be expressed linearly in 
terms of u, v y w, p, q, r. The forces are due to three sources: 1° the 
propeller thrust, 2° gravity, 3° the air. We shall assume that the 
propeller thrust (and moment, if any, arising from it) is constant; i. e., 
the motor is supposed to speed up or slow down under changed condi- 
tions so as to deliver a constant thrust. If 6 and <p are the small 
pitch and roll, the components of gravity are gO, —g<p, —g (see 
Bairstow, 144, 7u — w), and its moments are zero because the C. G. is 
taken as origin. The air forces and moments may be written as 
X, Y, Z, L, M, iVand developed as 

X = X 0 + X u u + X v v + X w w + X v p + X q q + X r r, (8) 

where X u , X v , are the " resistance derivatives" taken for 

the relative velocity of machine and wind. (X 0 and the propeller 
thrust cancel, so do Z 0 and g; Y 0 , L 0 , M 0 , N 0 vanish.) 

In the symmetric aeroplane half the resistance derivatives vanish 
and the six equations of motion separate into two sets of three each, 
one set for the longitudinal, the other for the transverse motion. 
These equations are (Bairstow, 148, 13 and 14 with 9 = 0) for longi- 
tudinal motion, 

du/dt =gQ + X u u 4- X w w + X q q, (9a) see ( 1 a) 

dw/dt =Uq + Z u u + Z w w + Z q q, (96) see ( 1 c) 

B/m. dq/dt= M u u-\ M w w+ M q q, (9c) see (26) 
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and, for transverse motion, 

dv/dt= -g<p- Ur+ Y 9 v+ Y p p+Y r r, (10a) see (16) 

A/m. dp/dt — E/m. dr/clt = L v v + L p p + L r r, (106) see (2a) 

C/m. dr/dt - £/m. dp/dt = A> + iV p p + N r r. (10c) see (2c) 

The integration of these equations gives the free oscillations of the 
aeroplane. 

Article 2. 
LONGITUDINAL MOTION IN SMALL GUSTS. 

A gust if not too severe may be treated by the method of forced 
oscillations. If the aeroplane is traveling on an irregular wind, we 
may regard the average wind velocity relative to the machine as that 
which should be used in the computation of the resistance deriva- 
tives, and we may regard the departures of the actual relative velocity 
from the mean as small quantities inducing additional forces into 
the equations of motion. 

Suppose first a head-on gustiness. This would introduce an extra 
term of the form X u u into the first equation, Z u u in the second, 
and so on. If, as a result of the 2ust, the machine tilted appreciably, 
the originally head-on gust would no longer be head-on, but would 
have components u u w x and give rise to the term X u u x + X w w x in the 
first equation. It is clear, nowever, that under the hypothesis of 
small oscillations, w x would remain small of the second order relative 
to i/ v The term X w w l could then be neglected relative to X u u lf 
unless X w much exceeded X u . 

We should in general allow a gust to have components u v v 1 ,w v p iy 
er n r, relative to the axes. This would take into account any possi- 
ble rotational motion in the gust. The rotational motion of a gust 
may be quite small. In the discussion by Glazebrook (Aeronautical 
Journal, July, 1914, pp. 272-301) nothing is accomplished relative 
to rotational gusts. I et it mav well be that the rotational element 
is of great importance. For the rotary derivatives, in the case of 
the machine whose derivatives are tabulated by Bairstow (loc. cit., 
159), are large. Thus a term M q q x = —210^! would be comparable 
with X u u x = — 0.14u t if g^were 1/700 of u^, i. e., if the gust were a uni- 
form whirl of radius 700 ieet. In the same way L p is large. In the 
machine that will be discussed in what follows M q is also large, 
viz., -150. 

The equations for the longitudinal motion in a general gust are 
(see 9a-c) 

du/dt -gd- X u u - X w w - X q q = X u u x + X w w, + X q q x . (11a) 

dw/dt - Uq- Z u u - Z w w - Z q q = Z u u x + Z w w x + Z q q t . (lib) 

B/m dq/dt - M u u - M w w - M q q - M u u x + M w w x + M q q t . ( 1 lc) 

The solution of these equations consists of two parts: 1° the so- 
called complementary function which gives the natural oscillations, 
2° the particular integral which gives the forced oscillations due to 
the gust. To effect a solution for the particular integral, we must 
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make some assumption as to the value of the components u lt w lt g t 
of the gusts as functions of the time. Before making such an assump- 
tion for the particular integral, the solution by the "operational" 
method may be indicated. (See Wilson, Advanced Calculus, p. 223.) 
Let D denote differentiation. The equations may be written 

(D-X u )u-X w w-(X q D + g)d =X u u, +X w iv 1 +X q q v (12a) 

- Z u u + (D- Z w )w - (Z q + U)D8 = Z u u x + Z a w, + Zgq v (126) 

- M u u- M„w + (k 2 B D 2 - M q D)9 = M u u, + M w u\ + M q q v (12c) 

where k\=B/m. These equations are solved algebraically by 
multiplying by the proper cofactor determinants and adding. Then 



D 



-X u -X w -(X q D + g) 
-Z u D-Z a -(Z q + U)D 
-M u -M,. 



M q D 



u = 



+ 



X u -X w -(X q D + g) 
Z U D-Z W -(Z q +U)D 
M u -M w k 2 s D 2 - MqD 

X w -X w -(X q D+g) 
Z W D-Z W -(Z q +U)D 
M w -M w l- 2 B D 2 -M q D 



w,(13) 



+ X„ 



' X u 



Za D-Z„, 



A 



M,„ k 2 s D 



(X q D + g) 
(Z q + U)D 



M q D 



or, if the determinant on the left be denoted by A, 



Au = 



+ D 



X u -X w - (XqD + g) 
Z tt D-Z w -(Zq+TJ)D 
JI4 -M w k 2 B D 2 -M q D 



M w 



- (XqD + g) 

k 2 B D 2 - MqD 



w. + 



X q — X w 

Zq D — Z w 

M„ - M,„ 



(14a) 



-9 
UD 



There are similar equations for w and 0, namely, 



Aw = 



D-X U X W - (XqD + g) 
-Z U Z W - (Z q +U)D 
-M U M W k\D 2 



M q D 



(146) 



+ D 



A0 = 



Zu ~(Zg+ U)D 

M u k 2 B D 2 — MqD 



u. + 



D-X u X q -g 
-Z U Z 3 - UD 
-M u M a k 2 B D 2 



2i» 



D — X u — X w 
— Z u D — Z w 



(14c) 



+ D 



D-Z w Z u 
~M W M u 



u x +D 



D~X U X w 



The general (literal) integration of these equations would be so 
complicated as to be useless. We shall make use of the formulas 
only after simplification by the insertion of numerical data. 
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Possible methods of treating gusts. — The only treatment of gusts 
which I have seen is that described somewhat popularly bv Glaze- 
brook (loc. cit.). He seems to state, as the main method of attack, 
that of small differences whereby it is assumed that the involved 
time over which the motion is to be studied is divided into small parts, 
and that the atmospheric conditions remain constant during each of 
these parts. By then regarding the differential equations of motion 
as equations in differences of the following form, 

Au' = (X'-wy + vV)A*, etc., 
A7^ = (mU + r[l\ 2 — q'h^At, etc., 

it is possible to compute, through a series of intervals At, the ap- 
proximate positions of the aeroplane. This method is, as Glazebrook 
states, exceedingly tedious, for At must be taken very small, indeed 
only a small part of a second in the case of a sharp gust, in order that 
the solution may be even approximately satisfactory for the differ- 
ential equations. Moreover, the whole calculation apparently has 
to be done from the beginning for each new type of gust whicn one 
desires to study. The method, however, is applicable in all gener- 
ality irrespective of the stability of the aeroplane. 

The reason that I have chosen to operate on the basis of small 
oscillations is that after a certain amount of preliminary calculation 
has been accomplished my formulas will enable me to treat very 
rapidly a series of very different types of gusts. My method is not 
applicable, of course, to machines which are not stable, for the oscilla- 
tions could not remain small with such machines, but it is probably 
doubtful whether the motion of the unstable aeroplane in a gusty 
wind is of very great importance, as the instability of the machine is 
not unlikely to cause indeterminately violent motions on relatively 
small gusts. I have tried to devise methods which would enable me 
to use graphical apparatus for obtaining the solutions here desired, 
but have been unable to throw the equations into a form which lends 
itself to such methods. 

Moreover, the coefficients which enter into the equations and into 
the solutions at all stages of the work are of such varying magnitudes 
that it is difficult to obtain any reasonably accurate results. It seems 
impossible — I have not yet succeeded in avoiding the difficulty — to 
eliminate the occasional necessity of subtracting numbers which are 
nearly equal in magnitude; thus the accuracy of the figures is, after 
subtraction, seriously impaired. As I was aware that the data fur- 
nished me were probably not accurate to three figures, I first made 
the calculations with slide-rule accuracy, only to find that the final 
results became wholly illusory, owing to the difficulty just mentioned. 
I have therefore had to recompute every thing with 4-place logarithm 
tables. Most of the figures which occur in the work are therefore 
4-place numbers. Those which appear to have only three significant 
figures generally have the fourth figure zero when occurring in 
formulas containing 4-place numbers. In the calculations toward the 
end of the research the 4-figure accuracy has become reduced to 
three or two figure accuracy, out it did not seem best systematically 
to reduce the numbers by the omission of two figures, although this 
reduction has occasionally been made in final calculations. 
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Article 3. 
NUMERICAL EQUATIONS FOR HIGH SPEED. 

The data for high speed are (see Hunsaker, p. 47): 
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X u = -0.128, X w = +0.162, X q = 0 

Z u = - 0.557, Z w = - 3.95, Z q = 0 

M u = 0, M w = + 1.74, M q = - 150 

5/m = fc 2 * = 34, 17= -115.5, £ = 32.17 

The cof actors 5 in the determinant A are — 

(Z q +U)D 



(15) 



D-Z w 

-Mu 



-M w 



lc 2 D D 2 -M q D 



7c 2 B D 2 — M q D 
-(X q D + g) 



Z> + 3.95 115.5X> 
-1.74 34D 2 +150D 



= 34Z> 3 + 2S4.3Z) 2 + 793.5P = 8 r 



- 1.74 34D 2 +150D 
-0.162 -32.17 



-X w -(XJD + g) 
D-Z w ~(Z q +TJ)D 



-(Z q +V)D -Z u 
Tc 2 D 2 — M q D -M u 



D-X u -(X q D + g) 
-M u l 2 B D 2 -M q 



~{X q D + q) D-X u 



* 5.508D 2 + 24.30Z> + 55.98 = 8 2 



-0.1G2 - 32.17 

D + 3.95 115.5D 
13A6D+ 127.1 = 6 31 

115.5Z? 0.5571 
34# 2 +150£> 0 
-18.94Z? 2 - 83.56^ = 5,0 



0.128 
0 



- 32.17 
34D 2 +150£> 



-(Z q +U)D 



— X u 

- M K 



D-X u 
-Z n 



D-Z w 



D-X u 



D-Z u 



34Z> 3 + 154.3J9 2 + 19.20D = 5 22 
- 32.17 D+ 0.128 

115.5D 0.557 
- 115.5D 2 - 14.782?- 17.92 = 5 32 



- 0.557 D+ 3.95 

0 - 1.74 

- 0.9692 = S 13 



-0.162 D + 0.128 
-1.74 . 0 
1. 742? +. 2227 -5 2 , 



LD + 0.128 -0.162 

0.557 D + 3.95 
# 2 + 4.078Z>+ .5957 = S 33 



The value of the determinant A is 

34Z> + 288.7D 3 + 833.0J9 2 + 1 15. ID + 3 1. 18 = 
34(Z> 4 + 8.490J9 3 + 24.502? 2 + 3 3S5D + 0.9 170). 
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(The value of the determinant checks by three calculations.) 
The roots of the equation 

f(D) =D* + 8A9D 3 + 24.5D 2 + 3.385D + 0.9 17 = 0 (16) 

determine the decrements and periods of the natural oscillations, 
and must be found. (Unfortunately these roots must be found with 
considerable accuracy, and the rough first approximations, such as 
are indicated by Bairstow, seem insufficient for our use.) Let it be 
assumed that one root is so large that it may be found approximately 
from 

D* + 8.49J9 3 + 24.5Z> 2 = D 2 + SA9D + 24.5 = 0. 

Then D= -4.245 + 2. 545i. 

If now r be an approximate solution of j(D)=0, a new approxi- 
mation may be had by assuming r + x, with x small, as a root. 
Then 

= /(**) r 4 + 8.497* + 24.5r + 3.385r + 0.917 
X jf'(r) 4T 3 + 25.47r> + 49r+ 3.385 

approximately. As r 2 + 8.49r + 24.5 = 0, the fraction simplifies to 

3.385r + 0.917 
X 23.08r + 211.4 • UW + - AU "> 

if r*» ~ 4.245 - 2.545i, This root of f(D) = 0 is therefore 

Z> = - 4.182 ±2.438i. 
The factor off(D) corresponding to this pair of roots is 

£ 2 + 8.3641? + 23.43. (17a) 

Let the other factor be D 2 + aD + b. Then 23.436 = 0.917 and 
b = .03914. Also, 8.364(.0391) + 23.43a =3.385 or 23.43a = 3.058 
and a = .1305. Hence the second factor is 

£ 2 + .1305Z> + . 03914. (176) 

As a check on the work we may multiply the two factors together; 
we find 

(D 2 + 8.364Z? + 23.43) (D 2 + . 1305Z? + .039 14) = 
D 4 + 8.494Z? 3 + 24.56D 2 + 3.385Z? + .9170. 

We can find, merely by careful trial, better factors as 

(Z> 2 + 8.359Z> + 23.37) (D 2 + .13082? + . 03924) = 

D* + 8.490Z? + 24.50J9 2 + 3.385Z> + .9170. ( 18 > 

The definitive roots of f(D) =A = 0 may therefore be taken as 

a= -4.180-2.430i, 6 = - 4.180 + 2.430i e ^ 

c= - .0654 -.1870*, a«-.0654 + .1870i ^19) 
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INTEGRATION FOR HIGH SPEED. 

The numerical equation for u is (see 14a) : 

34 (2? 4 + 8.49 2> 3 + 24.5 D 2 + 3.385 D + 0.917) u 
- (X tt 5 u + Z U S 21 ) fij + 2)5^ + Jf^j, 

= -34 (0.128 J? 3 +1. 160 D 2 + 3.385 2) + 0.917) (20a) 

+ 34 D (0.162 2? 2 + 0.715 D + 1.647) M> t 

-34 (59.37 D + 560.6) ft. 
The numerical equation for w is (see 146) : 

34 (2)* + 8.49 2) 3 + 24.5 2) 3 + 3.385 2) + 0.917) to 

= (XtA, + Z,^ + lfu,5 32 ) w 1 + 2)5^! + 1/^,2^! 

= -34 (3.95 2? 3 + 23.94 D*f 3.385 2? + 0.917) w t (206) 

-34 2> 2 (0.557 D + 2.458) u x 

+ 34 (509.5 2? 2 + 65.21 £ + 79.05) q,. 
The numerical equation for 6 is (see 14c) : 

34 (2? 4 + 8.49 2> s + 24.5 D 2 + 3.385 2> + 0.917) 0 

= Mfa& + T>d 13 u t + D8 2Z w l (20c) 

=-34 (4.412 2? 2 + 17.99 2) + 2.628) q, 

-34 (0.02851) 2^ + 34 2)(.05117 2> + .00655) w t . 
The solutions are of the type: 

u = C n #* + C 12 c* + C lz e« + <7 l4 c* + /«, 

w - C 21 c<* + tf 22 e<* + C^e* + C^c* + I m (21) 
0 - ^31^ + C Z2 e* + C^e* + C^e* + 2,, 

where a, 6, c, d are the roots of the biquadratic (see 19), Cy certain 
constants of integration, and I UJ I w , I 9 a set of particular solutions 
of the equations. We* shall determine I u , I w , I 9 in such a manner 
that they will not contain the functions e^, etc.; we may therefore 
determine in advance the relations between the twelve C's. (This 
will debar us from using as gusts u v w u q v those which are of the 
form Ce<*, etc.; but this restriction is not important — such a damped 
gust tuned to the damping and . period oi the machine is hignly 
improbable in nature.) 

If we substitute u, w y 6 in the equations (i4), the particular solu- 
tions must cancel out among themselves (since they can not cancel 
terms of the form e*) and leave 

(a - XJ C u eP - X w C n e* - (Xm + a) C zy e?" + similar terms = 0, 

-Z u O,c* + (a-ZJ 6 21 e«-(Z q + thaCL&*+„ =0, 

- MvC^ef*- M w C 2l eV + Qc* B D 2 -M q D) (7 31 ^+ =0. 
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These equations hold identically in t, and the coefficient of e^ et c. in 
each must vanish. The three homogeneous equations in the three 
unknowns C„, C 21 , 0 31 (or the simdar equations in t 12 , V ?2 , v 32 , v 13 , 
T 2 CL; C u , \, d u ) are consistent because a (or b c, d) is a root of 
the determinant A, and the solutions are: 



-X w -g 




— g a— X u 




a— X u — X w 


a-Z w - TJa 




— TJa —Z u 




— Z u a — Z w ' 



C n : C 2 i* O sl — 

with CU: G 22 . <7 M determined by the same functions of b. In .words: 
To obtain the ratios of the coefficients o ^ m u v, w, substitute 
X> = a in the determinants 8 n , 5 32 , 5 33 . Or G 12 : 6 21 . C 31 as 
13 46a + 127 1 : - 1 15.5a 3 -14.78a- 17.92 : a 2 + 4.078a + .5957 
ov T: <£! <? 31 = 13 46a + 127.1 : 950.8a + 2560 : -4.281a-22.81. 
This gives C n : C 2l : O n as 

70.8-32.7 *:- 1414-2310 i:- 4.92 + 10.40 i or as 

I- _4.04-34.52i: -. 1 132 + . 0946 i. 
The values of C l2 : C 22 : are th| conjugates 

! . _ 4.04 + 34.5i: -.1132 - .0946 i. 
To find C 1S : C 23 : (L we must substitute c = - .065 - .187 i in the same 
determinants. Then 

C ■ C: C„ = 13.46c + 127.1:. 33c -13.39 :3.947c + .5565. This gives 

^ 13 * 23 33 

C 13 : C 23 : c/jj as 

126.2 - 2.516 i: - 13.37 - .0623 i: .2983 - .7380 i 
or l : -.1058-. 002587 i: .002478- .005799 %. 

The values of the conjugates are: 

O u : C 2l : C 34 = 1: -.1058 + . 002587 i: .002478 + .005799 i. 

The general solutions of the equation of motion are: 

ott .' M (22a) 
w= ( 2! 4.04 -34.5 i)C u e*+(- 4.04 + 34.ot)tf ia £ 

+ ( - .1058 - .002587 i) <V + ( - -1058 + .002587 %) C u e* + 

0 = (-.1132 + . 0946 »)^+(-. 1132- .0046*)^ (22c) 
+ (.002478 - .005799 i) <7 13 e* + (.002478 + .005799 *) tf 14 «* + 

From these equations we see that the heavdy damped short period 
oscillation (roots a, b) is about 34* times as strong in ^ as in u; 
whereas the mildly damped long period oscdlation (roots c d). U 
about 94 times as effective in u as m w. Moreover, the short period 
motions in u and w are about quartered- but the long period mot-ons 
Ze in opposite phase. The amplitude of the short period mo a m 
6 is aboutUr that of w; hence for each foot-second of short oscdlation 
•mu> thereabout 1° in 9. The amphtude of the long period motion 
in 6 is about .006 of that in u; hence for each foot-second of long 
oscillation in u there is about J° in 6. The damping of the short 
oscdlation is so strong that the amplitude is reduced to about one- 
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rtieth in one second where in the case of the long oscillation 
the reduction is only to about nine-tenths of its original value in one 
second; the relative amplitudes in the cases of u, w, 6 are more 
important in the case of the long than in that of the short period 
oscillation because the latter is so quickly damped out that the 
swing; may not get well started. However, the extreme magnitude of 
the short period oscillation in w as compared with u indicates the 
possibility of relatively violent accelerations in w; indeed, it is the 
short period oscillation which may account for initial difficulties 
whereas the long period oscillation accounts for the progressive 
troubles, due to gusts. 

There remain to be determined the values of the constants 0 of 
integration from the initial conditions of uniform flight, i. e., u = w = 
0 = q = O. Let the particular solutions have the initial values 7^, 
I™, ho- Then 

0 = C u + C i2 + C l3 + CL 4 + Iya, 

0 = ( - 4.04 - 34.5i) <?„+(- 4.04 + 34.5 i) C 12 

+ ( - .1058 - .002587 i) C l3 + ( - .1058 + .002587 i) G u + l w0 , 
0 = ( - .1132 + .0946 i) C n + ( - .1132 - .0946 i) C l2 

+ (.002478 - .005799 i)C l3 + (.002478 + .005799 i)C u + I 0o . 
0 = ( - .1132+ .0946 i)aC n + ( - .1132- .0946 i)bC 12 

+ (.002478 -.005799 i)cC 13 + (.002478 + .005799 i)dC u + r d0 . 
or 0 = (.703 - .205 i) C u + ( .703 + .205 i) C l2 + ( - .001246 - .000084 i) C n 

+ ( - .001246 + .000084 i) C u + 7 V 

The values of C n , C 12 and (7 13 , C u are conjugate imaginaries; hence 
C U + C» = A, C l9 +C U = B, i(C l2 -C n ) = C, i(C u -C lz ) = D are real. 
The equations may therefore be written 

0 = A + B + I uo 

0= -4.04 A + 34.5 1058 B + .002587 D + I wo 
0= -.132 ^i-. 0946 C+ .002478 B+. 005799 D + I 9o 
0 = .703 A + .205 C- .001246 B + .000084 D + 1\ 0 . 

The values for A, B, C, D are (as found by determinants and checked 
by substitution): 

A= - .0008856 7^+. 008198 7^+. 01621 7^-1.372 
C= - .003196 I uo -. 02803 7^+. 01476 7, 0 -.1543 7'* 0 , (23) 
£ = -(1-. 0008856)7^-. 008198 7^-. 01621 7^+1.372 V 9o , 
7>=.35777 WO -.2940 l wo - 172.0 l 6o - 29.89 V Qo . 

The solutions (22) of the equations of motion of the aeroplane in- 
volve imaginary numbers from which they may be freed by using 
A,B, C,D in place of 0 1V C 12 , O w C u . The equations then become 

u = e -*.m cos 2.43* + 0 sin 2.43*] 

+ e-- 0G54t (B cos .187* + 7? sin .1870 + I u , 

w = e- 4A8t [(34.5 C-4.04 A) cos 2.43* 

-(34.5 A + 4.04 O) sin 2.43*] 
+ e~ MUt [(.002587 D- .1058 B) cos .187* 
-(.002587 B+. 1058 D) sin .187*] + I w , 
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0 = 6-<"«[-(.1132^ + .O946 C) cos 2.43* 

+ (.0946 JL-.1132 C) sin 2.43*] 
+ cr-°«« [(.00278 B + .005799 D) cos .187* 
+ (.002478 D- .005799 B) sin .187*] + /, 

These formulas enable us to study any particular gust we desire. 

It is merely necessary to find the particular solutions, then the 
constants A, B, O, D. We shall reduce the coefficients in the paren- 
theses. Then 

u==e -4.m cos 2.43*+ Q sin 2.43*) 

+ er-™* (B cos .187* + !? sin .187*) + I u , (24a) 

w = e- 4ASt (A' cos 2.43*+ C sin 2.43*) 

+ «r-°««* (B' cos .187* + 7)'sin .lS7*) + 4, (246) 

0==€ -4.i8t cos 2.43*+ <7" sin 2.43*) 

+ er-«»« (£" cos .187* + 7)" sin .187*) + 7*. (24c) 

where 

A' = - .1066 7,0-1.0001 7*0+. 4436 7* o +.220 7'* 0 , 

C' = . 04346 7^-. 1696 7^-. 6190 7* 0 + 47.93 I' eo , (25) 

# = .1066 7^ + . 000107 7„ 0 -.4436 7* 0 -.220 7'*, 

7?'= -.03523 7 tt „+.03112 7„ 0 + 18.20 7* 0 + 3.158 7'^ 

^1" = + .0004024 7^+ .001724 7 W0 - .003231 7* 0 + .1698 7V 
C"= + .0002778 7^-. 003947 7^-. 000136 7* 0 -.1123 7' 9o , (26) 
£" - - .0004024 I uo - .001724 7^- .99676 7^- .1698 7'*, ■ 
7?" = .006683 7^-. 000681 7^-. 4261 7, 0 -. 08201 7'* 0 . 

In any particular case the calculation of the coefficients in (24) 
from (23), (25), (26) is likely to be relatively simple because there 
are so many terms that for that case may be negligible. 

Article 5. 

SOME SPECIAL GUSTS. 

If we wish to represent a gust which, starting from the condition 
of still air, increases to a certain intensity J we may use the function 

J (I-*-*). (24) 

The value of r determines the sharpness of the gust. If r = l, the 
gust has reached about two- thirds of its value in one second; if r = 5, 
the gust has reached two-thirds of its value in one-fifth of a second; 
if r = the two-thirds intensity is reached in 5 seconds. We may per- 
haps regard r = 1 as giving a moderately sharp gust, r = 5 as giving a 
very sharp, and r=% as giving a tolerably mild gust. The Function 
(24) has trie advantage of being in such form that the determination 
of the particular integrals is easy. (See Wilson's Advanced Calculus.) 
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Case 1. Head-on gust — mild. u^J ( 1 — e~- 2t ). 

In equations (20) we let u x = J ( 1 — e-- 2t ), w 1 =g 1 = 0. Then 



7 tt = - J (1-.247 er*), luo = -.753 J, 

7„ = .082 J e- 2 ', 7^ = - .082 J, 

1 Q « _ .00495/ e--» 7'* 0 = - 0049e/, 

7'* = .00099 J er* } r Bo = .00099J. 



(N. B. — The total increase J of the wind occurs everywhere as a 
factor and may be omitted — the results then are for an increase of 
1 foot-second.) 

u = Je-- 0fS54i (.622 cos .187*+. 630 sin .1870-/(1 -,247<r*), 

w = j€- 4 - 18< (-.004 cos 2.43*+. 003 sin 2.430 - Je~ MhAt {. 078 cos .187< + 
.059 sin .187<) + .082J<r- 2 ', 

^ = ^-065^( 00495 cos .1872-.0031 sin .187t)- .00495 Je- 2t . 

It appears from these equations that the effect of a mild head-on 
gust of magnitude J is as follows: (1) The machine takes up an easy 
slowly damped oscillation in u of amplitude about 89 per cent of J; 
after the oscillation dies out the machine is making a speed J less rela- 
tive to the ground and hence the original speed relative to the wind. 
(2) There is a rapidlv damped oscillation in w of rather small magni- 
tude and a slowly damped one of about 10 per cent of J, the final 
condition being that of horizontal flight. (3) There is a slow oscilla- 
tion in pitch of about .0058 J radians or about .32 J°. If the mag- 
nitude J is great, the pitching becomes so marked that the approxi- 
mate method of solution can no longer be considered valid — a gust 
of 20 foot-seconds causing a pitch of some 6°. As the period is long 
(about one-half minute) the pilot should have ample time to correct 
the trouble before it produces serious consequences. * 

The result of a tail-on gust is the opposite of that of the head-on 
gust and therefore need not be treated separately. For the head-on 
gust J is negative; for a rear gust, positive. 

To calculate the stresses on the machine or operator caused by the 
gust we have merely to find the accelerations du/dt and dw/dt of which 
the first is (approximately) — 

du/dt = Je- M5 %08 cos .187*-. 16 sin A87t)- .05Je- 2t . 

This acceleration reaches a maximum of something of the order of 
J/10; and if J should be 20 foot-seconds, the acceleration would be 
only about 2, or 6 per cent of g — not a large amount. The accelera- 
tion dw/dt is likewise small. (N. B. — The initial accelerations du/dt 
and dw/dt should vanish, because the gust starts from zero. That 
the initial values are not exactly zero in the above formulas is due to 
the roughness of the final calculations for u and w.) 
The path of the machine varies from the horizontal by the amount 
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which accounts for the effect of the vertical velocity and of the climb- 
ing in the path. The result is (roughly) 

g-JjV^CS cos .187*-. 4 sin .I87t)dt-.5e--*dt, 
2 = J[e-.0654*( COS .1872 + 3 sin .1870 + 2.5e- 2 <- 3.5]. 

The motion is oscillatory approaching as a limit z= —3.5 J. The 
machine wall rise 70 feet when the gust is 20 foot-seconds head-on. 

Case 2. Up gust — mild. w l = J{l — e~- 2t ). 

7 u -.305 Jfr** 9 Ao = .305 J, 

4 = J(1-1.012*- 2 '), 4 o =-.012 J, 

le - .000737 Je~ 2t , I* = .000737 J, 

7* = - .000147 Jer4f, T eo = -.000147 J. 

u = Je- 06 ' oU (- .305 cos .1872-. 0108 sin .1870 + .305 

w = Je~ 4ASt (-. 02 cos 2.432 +.026 sin 2.432) + Je- 065 %032 cos .1872 + 
.002 sin .1872) + J(l-1.012e- 2 0, 

0 = Je- OQ5it (.0008 cos 1872+.0017 sin .1872) + .00074?- 2 *). 

The effect of the up gust is to set up a small long oscillation in u 
of magnitude about 0.3 J, a very small oscillation in w, and a long 
oscillation of intensity .0018 J radians or .11 J° in 0. The compar- 
ative effects on the velocity and angle in the case of head-on and up 
gusts show that the up gust is only about one-third as effective as the 
head-on gust. The accelerations in the case of the up gust are all 
small. 

To find the displacement in a vertical direction we integrate as 
before. 



0 = ^(a;+ 115.50)eZ2. 



It is scarcely necessary to trouble with the trigonometric terms 
partly because the motion is less pronounced than in Case 1, partly 
because there is here the secular term J t, which will carry the machine 
up with the gust and will be the chief effect after the lapse of a short 
time. 

A down gust is in every way the opposite of an up gust and need 
not be separately treated. 

Case 3 . Rotary gust — mild. q A = J(l — e~- 2t ) . 

I u = - J(610.6-475.5e- 2 <), I u0 = -135.1 J. 
I v; = J(S6.21-74.87e- 2 '), 7^ = 11.34 J. 

Ie= J(2.865 + .691e-- 2 0, 4> = 3.556 J. 
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Z*« - .138 Je - 2t , I l eo= - .138 J. 

I u = Je- 4ASt (Ai) cos 2.43*+. 1875 sin 2.43*) 
+ Je- M54t ( 134.7 cos .187*- 059 sin .187*) 
-J(610.6-475.5e -*)* 

J w = J 6 - 4 - 18 «(4.()l cos 2.43*- 16.82 sin 2.43*) 

+ Je- 0654 <(- 15.95 cos .187*4-70.08 sin .187*) 
4- J(8G.21-74.87e- 2 0, 

= Je-^(-.0698 cos 2.43*4-. 0223 sin 2.43*) 
4e/6- 0654 *(-3.487 cos .187*- 2.414 sin .187*) 
4*7(2.865 4.691 e - 2t ). 

The effect, of the rotary gust is a long oscillation in u (the short 
one is negligible) of magnitude about 670 J, a short oscillation in w 
of about 17 J and a long one of about 71 J, a long oscillation in 6 of 
about 4.1 J. The comparison with former cases may be made by 
supposing first that the oscillation in u may reach some 20 foot- 
seconds. Then J =1/33 =.03. The amplitude of the oscillation 
in 6 is then some 0.12 radians, which is an amount comparable with 
the 6° of Case 1. To get an idea of what J=.03 means, we may 
note that if a gust of 20 foot-seconds is due to a whirl of the air as 
a solid body with ^ = .03, the radius of the whirl is 660 feet. We 
may therefore say that the effect of a whirl of radius 660 generating 
velocity of 20 foot-seconds is of itself about equal to that of a head-on 
velocity of that amount. If, however, a machine ran into such a 
whirl, it would experience both the effect of the whirl and of the 
linear velocity generated by it and would be disturbed considerably 
more than if it had encountered a pure head-on gust. We may 
therefore say that if the head-on gust arises from a whirl of mate- 
rially less than 660-foot radius, the effect of the whirl is quite con- 
siderably larger than that due to a straight head-on gust of equal 
magnitude. 

The conditions after enough time has elapsed to allow the expo- 
nential term to become small is 

I u = -610.6 J. 4 = 86.2 J. I e =2Mo J. 
It is therefore seen that the machine takes up the head-on velocity, 
acquires a small upward velocity, and is inclined at an angle 2.865J 
radians to the horizontal, these effects being due exclusively to 
the rotary motion of the air. The path in space could be obtained 
by integration, but (like the effects previously mentioned) would 
not be the true path if the rotary motion were accompanied by 
horizontal or vertical linear gusts. It seems therefore scarcely 
worth while to find the path. 

The value that I attach to this theory of rotary gusts does not 
arise so much from the fact that such gusts seem nowhere to have 
been treated as from the revelation of the powerful effects of such 
gusts. When a machine is flying low it must expect to meet air 
which has been set in rotation by the friction of the wind against 
the ground, against buildings, or against trees. It seems certayi 
that very material angular velocities might be set up and that these 
might (owing to their short radius) induce only moderate linear 
gusts. In such cases, if they can arise as assumed, the machine 
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might behave very much worse than could be foreseen when nothing 
is known of rotary gusts. It is not unlikely, however, that rotary 
gusts would be very irregular themselves and that, before the 
machine could feel the full effects of one, the crust might have dis- 
appeared. In the same way rotation could be generated at the 
interface between dark and light regions of air — indeed any sharp 
relative motion of the air is likely to contain rotation. 

Case 4 . Head-on gust — moderate. u 1 = J(l — e- t ). 

I u = -J r (l + .09876e- < ), 7 tt0 = -1.09876 J y 

I w = .1307 Jer*, I wo = .1307 J, 

I e = - .00196 J<ET*i Ie 0 = - .00190 J, 

l'e= + .00196 Je-<, I'o 0 = + .00196 J. 

u = Je- 4 - 18 ^ - .000676 cos 2.43*- .000486 sin 2.43* 
+ Je-- 0654 '( 1.09944 cos .187*-. 1528 sin .1870 

-J (1 + .09876<r'), 

w=Je- AA8t (-. 01405 cos 2.43*+. 02528 sin 2.43*) 
+ e/e- 0654 <(-.1159 cos .187*+. 01493 sin .187*) 
+ .1307JV-', 

e^Je- 4 - 18 ^. 0001207 cos 2.43*- .00000895 sin 2.43*) 
+ Je-- 065 %001838 cos .187*-. 006755 sin .187*) 
- .00196 Je-K 

The short oscillation in u is negligible not only in regard to its 
magnitude but even as far as accelerations are concerned. Then 

du/dt = Je-' 0Q5it (-.l cos .187*+. 21 sin .187t) + .Mer*. 

This is at most about .25 J, or 5 foot-seconds 2 if J = 20. The short 
oscillation in w is considerably smaller than the long, but when the 
coefficients -4.18 and 2.43 are brought in by differentiating to find 
dw/dt, whereas - .0654 and .187 are brought in by the long oscilla- 
tion, it appears that the short oscillation is effective in determining 
the acceleration. Thus 

dw/dt = Je- 4 - m (A2 cos 2.43*-. 07 sin 2.43*) 
+ Je- 065 %01 cos .187) -.13 Jer*. 

The amount of this acceleration is at most about J/12, one-third that 
in u; the effect, however, is produced very quickly, in the first half 
second. 

In integrating to find the path in a vertical plane we may neglect 
the short oscillation, because in this case we divide by —4.18 and 
2.43, whereas for the long oscillation we divide by —.0654 and .187. 
Then 

z = j t (w+U5.50)dt 

= jj\e- 0(i54t (.106 cos .187*-. 765 sin .187*) - .095er t ]dt 
= J r 6 - °« 5 «(2.3 sin .187* + 3.5 cos 1.87*) + .095 Jo 
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The final condition is a rise of -3.6 J, an amount which agrees with 
that in the case of the mild gust (Case 1) in as far as the rough calcu- 
lation of that case permits us to judge. 

Case 5. Up gust — moderate. w 1 = J(l — ert) . 

I u = .0773 J<r*, I u o= -0773 J, 

lw= ~ J (1 - e-t), l w0 = .205 J, 

7 0 = - .003069 Je~\ I 6o = - .003009 J, 

r 9 = .0030(39 Je-\ V 0q = .003009 J. 

u = Je~ 4 ' 18t (- .002041 cos 2.43*-. 00651 sin 2.430 

-f J<r- 0654 '( = . 07466 cos .187*+. 4034 sin .1870 + .0773 Jer*, 

w = Je- 4ASt (-.2139 cos 2.43*+. 1174 sin 2.430 

+ Je- m54t (. 008943 cos .187* - .02337 sin .187Q- 1.205*-*), 

e=J e -4-m(, 0009 148 cos 2.43*+ .000487 sin 2.430 

+ J<r-- 0654 '( + .002154 cos .187*- .001432 sin .187*) - .003069 Je"'. 

The short oscillation is negligible in u as far as concerns u itself. 
In calculating the acceleration du/dt the short oscillation is not 
negligible relative to the long; but the acceleration is small any way. 
The effect of an up gust J on u is about one-third the effect of an 
equal head-on gust (see Case 2). 

The short oscillation is the main thing in w — its amplitude is about 
J/4, whereas the amplitude of the long oscillation is about J/40, or 
one-tenth as much. The acceleration dw/dt may therefore be cal- 
culated exclusively from the short oscillation; it is 

dw/dt = Je- 4ASt (l.2 cos 2.43*)- J (l-^O- 
This means values approximately as follows: 

* = 0, %j M, >2, Ma 
acc. = 0, - .35 J, - .6 J, - .7 J, - .6 J. 

If J should be 20 foot-seconds, the maximum acceleration would 
be about g/2, even a gust of 10 foot-seconds would produce an accel- 
eration of g/4. Such accelerations coming upon the pilot in one-half 
a second might considerably surprise and disturb him. An addition 
of 25 to 50 per cent in the apparent weight of the machine could 
hardly strain it to an appreciable extent in view of the large factor 
of safety used in the design. (N. B. — For an up gust J is negative. 
For a down gust the operator would lose 25 to 50 per cent of his 
weight.) 

The path of the machine in space is not of great importance in 
this case. The chief feature is the general drift of the machine with 
the current. 

Case 6. Rotary gust — moderate. q } = J (1 — e~ l ). 

As we know so little of the rotation in the atmosphere and as 
nothing particular of interest seems to be indicated for this case 
over and above what was found in Case 3, we shall not carry out the 
calculations. 
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Case 7. Head-on gust — sharp. u x = J (1— r -5 *). 



I u = -e/(l + .01S72 e~*) h I uo = -1.01872 J, 

4= - .05102 Jr-« 4 0 = - .05102 J, 

/<, = - .0008896 Je-*, I 0o - - .0008890 J, 

I'e = 004448 Jer**, I{ 0 = .004448 J. 



u = j e -*.i*t(_ >0 05632 cos 2.43M- .003980 sin 2.430; 
hJ^- 085 «( 1.02435 cos .187*-. 3294 sin .1870, 
-J(l + . 01872 e- bt ), 
w = Je- 4AHt (.ld03 cos 2.43*+. 1782 sin 2.430, 

hJe-' 0654t (-AQ93 cos .18W+.0322 sin .1870, 
-.05102 Je~ 5t , 
e=Je- 4A * t (.Q0Q26 cos 2.43*- .000984 sin 2.430, 
+ J<r*- W5 %000628 cos .187*- .006755 sin .1870- 
- .6008896 Je~ 5t . 

Here again the short oscillation in u is insignificant. The long 
oscillation as in Case 4 has an amplitude a little in excess of J. The 
acceleration du/dt is small of the order J/5. The reason that a sharp 
head gust does not give a large value to du/dt is probably because 
the gust can blow through the machine; the acceleration is therefore 
not large except at the loops of the slow oscillation. 

The short-period oscillation in w has now become stronger than 
the long oscillation and the acceleration dw/dt is mostly due to it 
and may be written 

dw/dt - J <r"*( - .25 cos 2.43*- 1.13 sin 2.430 + .25 Je~ 5t . 

The value of the acceleration never gets large because it is damped 
out before the sine term gets effective — perhaps —0.4 J would be 
about its maximum value. A sharp head-on gust is therefore about 
half as effective as a moderate up gust of the same intensity. Since 
up gusts are perhaps not likely to be as intense as head-on gusts, we 
might hazard a guess that sharp head-on gusts would inconvenience 
the pilot about as much as moderate up gusts. 

The most important terms in the path in space are 

z = Je- 065it (l.2 sin .187* + 3.5 cos .1870-3.5 J. 
The total rise is again —3.5 

Case 8. Up gust — sharp. u\ = J( 1 — er H ) . 

I u = .06621 Je-*, I uo = .06621 J, 

4= -J( 1^.5605 6" 5< ), I wo = - .4395 J, 



u = Je- 41st (--057U cos 2.43*+ .006 sin 2.430 

+ </6- 0654 '(-. 00907 cos .187*+. 3285 sin .1870 
+ .06621 Je~ 5t , 

w = Jer*- l *(A378 cos 2.43*+ 1.947 sin 2.430 

+ Je- 065 %00181 cos .187*-. 03474 sin .187*) 
- J(l-.5605 e~ 5t ), 



f 9 = - .00778 Je~ 5 S 
'* o =.0389 Je~ 5t , 



I Bo = - .00778 J, 
0389 J. 
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0 = Je-- 4m (.OO59 cos 2. 43*-. 0122 sin 2.430 

+ Je- 0654 '(.001883 cos .187*+ .0008607 sin .187*) 
- .00778 Je~ 5t . 

The oscillation in u is of long period; and the acceleration in u is 
small. The oscillation in w has a short-period term of great impor- 
tance at the start, but except for this there is very little oscillation 
in u\ The acceleration is 

dw/dt = Je- 4A8t (2.9 cos 2.43*- 9.2 sin 2.43*) -2.8 Je~ 5t . 

(N. B.— The value of dw/dt when * = 0 should be 0 instead of J/10. 
The failure to check seems due to multiplication of errors, which is 
unavoidable. The accuracy of the work in Case 8 and Case 5 appears 
reduced to two figures.) The acceleration is now very serious indeed ; 
it is about -9.2 Je~ 4A8t sin 2.43*, as the other two terms come 
near canceling. The maximum value occurs when * = .217, a little 
over one-fifth of a second, as is then about —1.85 J. If J should 
be as large as —18 foot-seconds, the acceleration would equal gr?=35. 
Clearly such a sharp gust if it existed would be very dangerous from 
the sudden forces it would bring into play. As the machine, how- 
ever, would travel only about 24 feet during one-fifth second, it is 
reasonable to doubt whether in so short a distance so large a change 
in vertical air velocity could occur. 

The path in space is found to be approximately 

z= -1.2 Je- iA8t cos 2.43*+ 1.1 JerW* cos .187*-.! Je' :,t + .2 J-Jt. 

The final effect is the general drift with the gust, less a lag of J/5. 

Article G. 

THE CONSTRAINED AEROPLANE. 

If an aeroplane is constrained to remain always horizontal by 
mechanism which does not otherwise alter the machine or its dynam- 
ical properties, the equations of motion in a gust may be found from 
our previous equations by setting 0 = ^ = 0. Then 

(D - X u ) u - X w w = X u u x + X w w t + X q q ly 

-z u u+(D- z ll} ) w = z u u x + z w v> x + z,//, ; 

- M u u - M w w = M u u x + M w u\ + Mqq, + F y 

where Fis the effective force due to the constraint and is assumed to 
affect moments only, not component^ of horizontal or vertical force. 
The last equation merely determines F. 

With the numerical data we find for high speed 

(Z>. + 128W-.162w = -.128^ + .162^, 
.557u + (D + 3.95)w - - .557^ - 3.95w„ 

F= - .1740 + ^) + 150&. 
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The natural motion of the machine when slightly disturbed in 
steady air is found from 

k = \ D + ]S d + 3, = 1)1 + 4 - 078P + - 59s = °- 

The roots are 

• D= -2.039 ±1.887= -3.926 or -0.152. 

We thus find the first result: The machine, when disturbed, does 
not execute a double damped oscillation, but has an aperiodic motion 
of the form 

C e s™+C 2 e-°- 15t . 

The two damping factors -3.93 and -0.15 lie between the values 
— 4.18 and —.0654 previously found. 

The unconstrained machine was stable for the speeds 79, 51, and 47 
mile-hours; unstable for 45.2 mile-hours and lower speeds. If we 
take the data for 47 mile-hours and use them for the constrained 
mo1 ion, we find 



A"=D 



of which the roots are -1.51 and -I-.10. The natural motion of 
the machine is therefore of the form 

The second factor indicates instability; the motion due to it increases 
instead of subsides and reaches 2.78 times its original value in 10 
seconds. We thus find the second result: The machine, when con- 
strained, becomes unstable at a higher speed than when free — it is 
to this extent a more dangerous machine. 

We shall now return to the case of high speed and compute the 
effect of certain gusts on the constrained machine for comparison 
with the effect of the same gusts on the free machine. The general 
solutions are 

u = - .0426 C\ e~ 3 -^ + C 2 e-^ + I u , 

W =C X e- 3 - 93 '-.147 C 2 e~ A5t + I w . 

a =-.148 l uo - 1.006 I w0 , 

C 2 =- 1.006 / VO -.0429 I wo . 
A'u= -(.128 D+.598) ^+.102 Dw v 
A'w = -(3.95 D + .598) w t - .557 Du x . 

Case!. Head-on gust — mild. u x = J (\—e~- 2t ). 

I u =-J (i-f-3.20 e-*)i I uo = -4.20 J, 
4 =.622 Jc-v, / WO = .022 J. 

<a = 4.19 Je-'^.-J (1 + 3.19 «-*), 
w= - .62 Je-™ + .Q2 Je- 2t . 
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The machine takes up the gust as before, of course. There is no 
oscillation. There is practically no acceleration in either u or w. 
The path in space is 

z = J (4.i e - 15t -SA e~- 2t )-J. 

The total rise is only — J. In every way the motion in this case is 
easier in the constrained than in the free aeroplane. 
Case 2. Up gust — mild. w l = J(l—e — ' 2t ). 

£- -.186 j _ mn * 
4=-«7(l- 1.079 e -*), > Z org / ' 

W= - .052 J(?- 3 - 93 '-.027 Je-' 15 *- J(l- 1.079 €"•*). 

The motion is again exceedingly moderate in all respects. 
Case 3. Rotary gusts. These can have no effect except upon 
the constraining moment F. 

Ca se 4. Head-on gust — moderate . u l = J(l—e~~ t ). 

I u = - J(l + .1895 e~*). I uo = - .1895 J, 
7 U , = .2246 Je~\ I wo =.2246 J. 

u=.002 J>- 3 - 93 <+ 1.187 Je"- l5 < - J(l -f .189 «-«), 
w = - .05 Je~ 393t -.174: Je~ A5t +.224: Je~*- 
du/dt= - .008 Je~ 3 -* 3t -.180 Jtf-- 15 <+1.89 
dw/dt = . 197 Je~ 3 - 93i + .027 Je-- 15 '-.224 
2=1.16 Je~ Abt -.22 Je-t-,94: J. 

The motion is again decidedly moderate. 
Case 5. £7p — moderate. w x = «/(l — e~ l ). 

I u = - .0053 J*-', 7^= - .0653 J, 

/„ = - J( 1- 1. 350 e~ l ), /«,«>=. 350 J. 
^=.0144 J r ^- 3 -»« + .0507 e- - 15 '-.0653 Je-', 
w = - .343 Je- 3 - 93 '-.007 e— 15 *~ J(l-L350 
<2w/(ft = + U35 Je~™ 3 - 1.35 

The motion is easy except for the acceleration in which has a 
maximum when t = .46 and is then equal to about — .62 J. If the 
gust should have an intensity of 10 foot-seconds the maximum 
acceleration would be about g/5. 

Case 6. Head-on gust — sharp. u l = J(l—e~- st ). 

I u =-J(l+ .00795* " 5 ), 4 o =-1.008J, 
/„, = - .5275 Je~' ot , I w0 = - .5275 J. 

u=- .029 Je ~ 3 - 93 * + 1 .037 Je - J ( 1 + .008 e - 5t ) . 
io = . 680 Je" 3 - 93 '-.152 Je"- 15f -.528 Je~ 5t . 
dw/dt= -2.67 Je~ 3 - 93 ' + .02 Je—w+ZM Je~ 5t . 
z = -.173 Je- 3 - 93 < + Je-- 15 < + .103 Je~ 5 <-.93 J. 
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The motion, including acceleration, is moderate. 

Case 7 . Up gust — sharp. w l = J(\—e~ ht ). 

/„ = .153 Je~ 5t , I U0 = A53 J j 



I w = -J(l + 3.628e- 5 <), I wo = -4.628 J, 
u = - .197 Je-*™ + .04:4: Je- A5t + .U,Z fer-fy 
to « 4.634 Je-*->*-.006 +3.628 e~ 5t ). 

dw/dt= -18.2 J<*- 3 -™ + 18.2 Je" 5 ', 

2= -1.18 j6- 3 - 93 < + .04 Je- - 15 < + .73 ^- 5 < + .41 J- J*. 

The acceleration dw/dt has a maximum when £ = 5/11 when it is 
1.44 J. This is somewhat serious if J is 10 foot-seconds. 

We may now calculate roughly the moment F necessary to pro- 
duce the constraint. 



The last term is effective only when the machine encounters rotat- 



ing air and will be neglected here. 

Case 1. F=.\\ J{e—' ht -e~' 2t ). 

Case 2. F=J(.009 *- 393 < + .005 e-- 15 <-.014 e~- 2t ). 

Case 4. F=J(.009 <?- 3 - 93 ' + .030 e- ^-.039 e~t). 

Case 5. F=J(.06 e~*** + >0 012 e— 1S *-.0612 

Case 6. F= J(-.119 *- 3 - 93 < + .0266 ^- J >< + .0924 e~v). 
Case 7. J(-e-- m + e~ 5t ). 



I have indicated the general method, based on the theory of small 
oscillations, whereby the equations of motion of a stable aeroplane, 
whether free or constrainecl to fly without pitch, whether in steady 
or gusty air, may be completely integrated in such form that, after 
a certain amount of preliminary calculation, the effects upon the 
motion of a large number of different gusts may be determined with 
relative ease. So far as I am aware, no actual method of integrat ion 
nor any quantitative results of such an integration has previously been 
published with the exception of the descriptive popular lecture of 
Glazebrook cited above. I have carried through the actual deter- 
mination of the effects of gusts in the following cases: 

Head-on gusts rising from 0 to J feet per second with various de- 
grees of sharpness. 



Rotary gusts of the same type. 
Rear gusts and down gusts are included by merely changing the 
sign of J. For convenience, it has been assumed that the machine 
is in still air except for the gustiness; as a matter of fact gusts are 
usually superposed upon a general steady wind of other than zero 
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average velocity; but the conditions of flight in still air and in steady 
air are nearly identical, the only difference being that in the equa- 
tions of motion the resistance derivatives are calculated from the 
relative wind, whereas U is the actual velocity over the ground. 

It has been found that a stable machine, with controls untouched, 
running into a head gust of various sharpness and of total intensity 
J foot-second will swoop up, with some oscillation of no serious char- 
acter, to a new level about 3.5 J feet higher than its previous level. 
The constrained inachine will rise without oscillation to a new level 
only J feet, or a trifle less, higher than before. The path in a ver- 
tical plane is indicated in the diagrams drawn for me by Mr. T. H. 
Huff. The accelerations arising m the motion are not serious for 
either the machine or the pilot. It has been found further that a 
rotary gust may have considerable effect — though in the absence of 
data as to the intensity and regularity of rotation in the air no definite 
results can be formulated. Furthermore we find that up gusts 
operate chiefly in lifting the machine, whether free or constrained, 
with the gust. The path in space is given in the diagram. There is 
here in the case of sharp gusts a considerable momentary acceleration 
in the vertical which may reach a magnitude of about 1.5 J foot- 
seconds. 2 This would not seriously stress the machine, which is 
designed to stand accelerations of 6 g to 8 g in maneuvering, but 
owing to its sudden and unexpected appearance tins acceleration 
might incommode the pilot — it is indeed trie familiar phenomenon of a 
"bump." 

It follows, therefore, that the introduction of the constraint, 
whether by gyroscopic or other means, serves only to eliminate the 
natural oscillation in pitch and to diminish, in the case of the head 
or rear gusts only, the final change of level. As a rear gust of 20 
foot-seconds is found to drop the uncontrolled machine by more than 
80 feet in 15 seconds, flight at low altitudes is more dangerous in the 
unconstrained than in the constrained machine. However, the 
elapsed time is sufficiently great to enable the pilot to check the dip 
by a suitable movement of his elevator. 

To offset any advantages derived from the constraint, we find that 
this particular machine, when constrained, becomes unstable at a 
speed between 47 and 51 mile-hours, whereas the free machine remains 
stable down to a speed between and 47 mile-hours. 

Massachusetts Institute of Technology, 
Boston, Mass., October 7 , 1915. 
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